Global study of the spectroscopic properties of the first 2+ state in even-even nuclei 
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We discuss the systematics of the 2"*" excitation energy and the transition probabihty from this 2"*" 
to the ground state for most of the even-even nuclei, from ^''O up to the actinides, for which data are 
available. To that aim we calculate their correlated J = Q ground state and J = 2 first excited state 
by means of the angular-momentum and particle-number projected generator coordinate method, 
using the axial mass quadrupole moment as the generator coordinate and self-consistent mean-field 
states only restricted by axial, parity, and time-reversal symmetries. The calculation, which is 
an extension of a previous systematic calculation of correlations in the ground state, is performed 
within the framework of a non-relativistic self-consistent mean-field model using the same Skyrme 
interaction SLy4 and a density-dependent pairing force to generate the mean-field configurations 
and mix them. To separate the effects due to angular-momentum projection and those related to 
configuration mixing, the comparison with the experimental data is performed for the full calcula- 
tion and also by taking a single configuration for each angular momentum, chosen to minimize the 
projected energy. The theoretical energies span more than 2 orders of magnitude, ranging below 
100 keV in deformed actinide nuclei to a few MeV in doubly-magic nuclei. Both approaches sys- 
tematically overpredict the experiment excitation energy, by an average factor of about 1.5. The 
dispersion around the average is significantly better in the configuration mixing approach compared 
to the single-configuration results, showing the improvement brought by the inclusion of a dispersion 
on the quadrupole moment in the collective wave function. Both methods do much better for the 
quadrupole properties; using the configuration mixing approach the mean error on the experimental 
B{E2) values is only 50%. We discuss possible improvements of the theory that could be made by 
introducing other constraining fields. 



I. INTRODUCTION 

Self-consistent mean-field methods (SCMF) are the 
only computationally tractable methods which can be 
applied to medium and heavy nuclei and have a well- 
justified foundation in many-body theory [||. Recently 
there has been considerable progresses in using these 
methods to compute nuclear mass tables 0. One of 
the appealing features of the SCMF is that the proper- 
ties of all nuclei are derived from a fixed energy density 
functional that depends on a small number of universal 
parameters, and can be used for the entire chart of nu- 
clei. The Skyrme family of functionals which is used in 
the present study, depends on about 10 parameters for 
the particle-hole interaction with 2-3 extra parameters 
for the pairing interaction. 

In a previous study [1, [1] j we have used an extended 
SCMF theory to calculate the ground state binding en- 
ergies of the about 600 even-even nuclei whose masses 
are known experimentally. In particular, two exten- 
sions were introduced to treat correlation effects going 
beyond a mean field approach. The first is a projec- 
tion of the SCMF wave functions to restore symmetries 
broken by the mean field: particle numbers and angu- 
lar momentum. The second is a mixing of projected 
mean-field states corresponding to different intrinsic axial 
quadrupole deformations. These calculations were per- 
formed with the same energy functional as for the de- 



termination of the mean- field configurations, so they do 
not require to introduce new parameters. Our main aim 
was to determine the effect of correlations on masses. 
In particular, the error on experimental masses in mi- 
croscopically based methods presents arches between the 
magic numbers. The correlations added in Refs. [3, Q 
clearly reduce the amplitude of the arches in the mass 
residuals, but do not remove them completely. For the 
parameterization we have used, however, the arches are 
much more pronounced along isotopic chains than along 
isotonic chains, which suggests that their appearance is 
not only related to missing correlations, but also to defi- 
ciencies of the currently used effective interactions. There 
is a clear improvement when looking at mass differences 
between neighboring nuclei around magic ones, in partic- 
ular when crossing proton shells. Similarly, the system- 
atics of charge radii is also improved, particularly in the 
transitional region between spherical and well-deformed 
nuclei. Altogether, this study clearly confirms the im- 
portance of incorporating some beyond mean-field corre- 
lations explicitly in the method and not heuristically in 
the energy density functional. 

In this work, we extend our previous study to two new 
observables: the excitation energy of the first 2+ state 
and the B{E2) value for the transition between the first 
2+ and the ground state. 

A simple approximation that can be systematically ap- 
plied is to start from a set of constrained SCMF config- 
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urations corresponding to different axial quadrupole mo- 
ments and to project them on angular momentum. Then, 
one searches the projected configuration that has the low- 
est energy for each angular momentum. We will call 
this procedure minimization after projection (MAP). A 
more sophisticated procedure requires an additional step: 
for each J-value, the total energy is further minimized 
by mixing projected SCMF configurations correspond- 
ing to different deformations. The mixing of constrained 
SCMF configurations is called the Generator Coordinate 
Method (GCM), and is based on the solution of the Hill- 
Wheeler (HW) equation. We shall label configuration 
mixing calculations by HW.. Although numerically de- 
manding, this approach has nowadays been used to study 
the detailed structure of low-energy collective excitation 
spectra of nuclei up to the actinide region using non- 
relativistic Skyrme [1, i, 0, H, H and Gogny [l^Jlll, Jll 
interactions as well as relativistic Lagrangians [13|, Il4| . 
Here, however, as in [1, we aim at something differ- 
ent: the calculation of a few very specific properties of 
the lowest 0+ and 2+ states for several hundred nuclei. 
For this purpose, the numerical procedure necessary for a 
detailed study by the GCM is too costly to be applied on 
such a large scale with present computational resources. 
Our bias on lowest collective 0+ and 2+ states, however, 
permits to set-up a tailor-made numerical scheme that 
reduces the effort considerably. For the angular momen- 
tum projection, we will generalize the topological Gaus- 
sian overlap approximation [l^, [l^ (GOA) used in our 
previous work. The GCM calculations will also be re- 
duced in size by truncations of the configuration space. 



time-reversal and having the same occupation number, 
which limits the description to even-even nuclei, and non- 
rotating states. The only other restriction on the wave 
function is that the Slater determinant of the orbits is 
invariant with respect to parity and axial rotations. For 
this work we take the SLy4 Skyrme parameterization [l9[ , 
the same as we used in our previous global study. For 
the pairing interaction we choose an energy functional 
that corresponds to a density-dependent zero-range pair- 
ing force, with cutoffs at 5 MeV above and below the 
Fermi energy, as described in [20| . As in earlier projected 
GCM studies, the pairing strength is taken to be —1000 
MeV fm^ for both protons and neutrons. 

To avoid a breakdown of pairing correlations for small 
level densities around the Fermi surface, we enforce the 
presence of pairing correlations using the Lipkin-Nogami 
(LN) prescription as described in Ref. [21]. However, we 
emphasize that the LN prescription is only used to gener- 
ate wave functions of the BCS form; physical properties 
are calculated with the code promesse [l^], which per- 
forms projections on proton and neutron particle num- 
bers and provides the matrix elements needed for angular 
momentum projection. 

Mean-field states with different mass quadrupole mo- 
ments are generated by adding a constraint to the mean- 
field equations to force the intrinsic axial quadrupole mo- 
ment q to have a specific value. Higher-order even ax- 
ial multipole moments are automatically optimized for 
a given mass quadrupole moment. A typical calculation 
for a nucleus involves the construction of about 20 SCMF 
wave functions that span a range of deformations suffi- 
cient to describe the ground state. 



II. CALCULATIONAL DETAILS 

As mentioned in the introduction, calculations are per- 
formed along the lines of Ref. Q. We will briefiy sum- 
marize the most important points and give details only 
for the necessary extensions to calculate 2+ states and 
matrix elements of the quadrupole operator. The calcu- 
lations reported here go beyond mean field in three re- 
spects: (i) projections on good particle numbers; (ii) pro- 
jection on angular momentum J = and 2; and (iii) mix- 
ing of states with different intrinsic deformation. All the 
results presented in this paper include particle-number 
projection and we drop explicit reference to particle num- 
ber from the notation. 



A. DFT Calculations 

We use the code ev8 (see Refs. [13, [lH) to solve 
the SCMF equations for an energy functional based on 
the Skyrme interaction. The single-particle orbitals are 
discretized on a three-dimensional Lagrange mesh cor- 
responding to a cubic box in coordinate space. The 
code imposes time-reversal symmetry on the many-body 
state, assuming pairs of conjugated states linked by 



B. Beyond mean field 

Formally, eigenstates | JMq) of the angular momentum 
operators and Jz with eigenvalues J{J+1) and M are 
obtained by application of the operator 

2 7+1 f^^ f^^ 
Pmk^^^]^ daj^dp smiP)J^ d-f V*Jk R, (1) 

on the SCMF states \q). The rotation operator R and the 
Wigner function "Dfjj^ both depend on the Euler angles 
a,/3,7. 

In a further step, we consider the variational configu- 
ration mixing in the framework of the Generator Coor- 
dinate Method. Starting from the ansatz 

iJMfc) =^/jfc(g)|JM(/), (2) 
g 

for the superposition of projected SCMF states, where 
k labels the states for given J and M, the variation of 
the energy {JMk\H\JMk) /{JMk\JMk) leads to the dis- 
cretized Hill- Wheeler- Griffin equation 

[Hj{q, q') - Ej,k I.,{q, q')] fjMq') = « (3) 
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that determines the weights fjk{q) of the SCMF states 
in the projected collective states, and the energy Ej ^ of 
the collective states. 

We have to calculate diagonal and off-diagonal matrix 
elements of the norm and Hamiltonian kernels. For the 
sake of simple notations, we use a Hamiltonian operator 
in all formal expressions, although there is no Hamil- 
tonian corresponding to an effective energy functional. 



In practice, as it is common procedure we replace 
the local densities and currents entering the mean-field 
energy functional with the corresponding transition den- 
sities. 

The axial symmetry of the mean-field states allows to 
simplify the 3-dimensional integral over Euler angles to 
a one-dimensional integral: 

Ijiq,q') = {JMq\JMq') 



(4) 



Jyjqjyjq' Jo 



Hj{q,q) = {JMq\H\JMq') 



AfjqKjq' Jo 

with normalization factors 

»i 



dcos(/3)di(o(/3) {q\RpH\q') 



(5) 



Af% = (2J+ 1) / dcos(/3) dUf3) {q\Rp\q)- (6) 



For the calculation of B{E2) values and spectroscopic 
quadrupole moments we also have to evaluate matrix el- 
ements of the quadrupole operator, which is outlined in 
appendix [X] A more detailed discussion of the method 
can be found in Refs. [1, Q and references given therein. 



1. topGOA overlaps and Hamiltonian matrix elements 

In Ref. [l^l, we found that for the description of the 
properties of the ground state the J — projected over- 
laps can be computed with a sufficient accuracy for our 
purpose with a 2- or 3-point approximation to the inte- 
gral using an extension of the gaussian overlap approxi- 
mation called the topGOA J24|] . There the rotated over- 
laps are parameterized by 

r (g|g')e-'=(9'9')-P('3) 
{q\RpW)t = <^ or (7) 

where F(/3) — sin^(/?) and the subscript t specifies the 
topGOA approximation. This form satisfies the require- 
ment of the GO A that F{P) {3'^ for small (3 as well 
as the topological requirement that F{-k — fi) — F{(3). 
Projected matrix elements of the Hamiltonian are also 
needed; these were calculated assuming the functional 
form 







{q\RpH\q')t2 = {q\q')e-^'^^-^"^^'^^\h^{q,q')-h2{q,q')Fm 

{q\RpH\q')t3 = {q\q') e-^(«'«')^(«-'^(«'«')^'('^) [/lo(g, <z') " h,{q, q')F{(3) ~ h^{q, q')F\p)] 

I 



(8) 



for the 2- and the 3-point approximations respectively. 
In general, the 2-point approximation is adequate for 
heavy nuclei and large deformations, but the 3-point ap- 
proximation is necessary to describe light nuclei. We 
take points at [3 equal to zero, and to a value where 
{q\Rl3\q') ~ 0.5 for the 2-point approximation. A third 
point is added at /3 = 7r/2 for the 3-point approximation. 
This is important for matrix elements between oblate and 
prolate configurations, which have their maximum value 
at 7r/2, 

In Fig.m we show an example of an energy curve deter- 
mined by this procedure for ■^^Ar. Angular momentum 
projection changes the quadrupole moment from the in- 
trinsic one to the one observable in the laboratory frame, 
which now depends on angular momentum. Most no- 
tably, it is zero for J = states independent of the 
deformation of the intrinsic configuration. As a conse- 
quence, it is more convenient and intuitive to use the 
intrinsic quadrupole moment of the SCMF states to la- 



bel the projected states. The marked points correspond 
to the q values of the SCMF configurations that were 
previously calculated. They are connected with lines to 
distinguish the J = and J = 2 curves. 

The J — Q curve has two very shallow minima at 
deformations q w -100 and -hlOO fm^. The J = 2 
curve has a pronounced oblate minimum at —125 fm^ 
and a shallow secondary minimum at -t-175 fm^. For the 
MAP calculation, we next estimate the quadrupole mo- 
ment at the minimum by interpolating between the cal- 
culated points. We then redo the calculations at the esti- 
mated minimum to find the MAP energy and quadrupole 
properties. For ^^Ar, the minimum for J = is at 
qa = -96 fm^ with and energy Eq = -332.32 MeV. The 
corresponding quantities for J = 2 are (72 = —120 fm^ 
and E2 = -328.33 MeV. The MAP excitation energy 
is the difference, E20 = E2 - Eq ^ 3.9 MeV. This is 
80 % higher than the experimental excitation energy of 
2.17 MeV. This is a rather extreme case, in that the 2+ 
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FIG. 1: Energy landscapes for J = and J = 2 angular- 
momentum projected states in ^^Ar. The open circles show 
the q values of the calculated configurations, with lines drawn 
to guide the eye. Solid circles are ones included in the mixed 
configuration calculations for the global survey. 



of '^*Ar is very likely better described as a broken-pair 
two-quasiparticle state than as a field-induced deformed 
state. We will return to this point later. 



C. Matrix elements of the quadrupole operator 

The calculation of the quadrupole moments of pro- 
jected states requires the calculation of all components 
of the quadrupole tensor. Q2±i and Q2±2 are of course 
exactly zero for axial mean-field states with the z axis 
as symmetry axis as chosen here, but they have non- 
vanishing transition matrix elements between a rotated 
and an unrotated state. 

The detailed expressions for the quadrupole operator 
and its projected matrix elements can be found in ap- 
pendix]^ For axial states, as used here, only the matrix 
elements of Q20 and the real part of the matrix elements 
of Q21 and Q22 need to be calculated, which simplifies 
the computational task. 

To calculate matrix elements of the quadrupole oper- 
ator Q2fn some modifications of the GO A parameteriza- 
tion are necessary since the functional behavior of the 
operator depends on its azimuthal angular momentum 
fi. In particular, for the matrix element of Q2±i, the 
form F{(3) — sin^(/3) used in the polynomial expansion 
of Eq. ([8]) is not topologically correct. We therefore de- 
fine a topGOA by taking (i^o(/^) ^'^^ argument of the 
polynomial expansion 



There is an additional complication compared to the 
norm and Hamiltonian kernels. While for these scalar 
operators the kernels ([7]) and ([5]) are invariant under 
exchange of \q) and \q'), this is not the case for the 
quadrupole operator, where (?|-R/3Q2/j|<z') is not equal to 
{q'\Rl3Q2fj.\q)- To avoid the exphcit calculation of both, 
we express the latter matrix elements as a weighted sum 
of the former using angular-momentum algebra and the 
symmetry properties of the SCMF states A separate 
topGOA is then set up to calculate the projected ma- 
trix element from the (?'|-R/3Q2/^|'?)- It has to be noted, 
however, that the difference plays a role only for transi- 
tion matrix elements between states with different angu- 
lar momentum. As we are interested here in 2+ — > 0^ 
transitions only, the topGOA for matrix elements with 
exchanged arguments is needed for ^ = only, Eq. (IA5|) . 

An example for the fits of the integrand is shown in 
Fig. Ofor ^^Cr at a deformation of q — q' = 150 fm-^. 
Starting from the bottom panel, the three panels show 
the rotated overlap matrix element {q\Rf3Q2fj.\q) for /i = 
0, 1 and 2 respectively. The open circles are the points 
used to evaluate the integral by a 12-point Gaussian 
quadrature, as it would be used in a calculation for the 
complete low-energy spectroscopy of this nucleus. The 
solid circles are the points used for the topGOA fit, the 
resulting curves being indicated by lines. The agreement 
is excellent; the error associated with the topGOA is typ- 
ically less than 1 % for the /i = matrix element. This 
is the only one needed to calculate the B{E2) value of 
the 2+ O"*" transition (see appendix \K\ . The middle 
panel shows the integrand for the Q21 operator. In effect, 
only the middle point can be used for the fit because the 
integrand vanishes at /3 = and tt. Nevertheless, this 
approximation works rather well. It is less accurate for 
some non-diagonal matrix elements, particularly for ma- 
trix elements connecting configurations with very differ- 
ent deformations which are needed to describe soft nuclei. 
The top panel shows the matrix element for the operator 
Q22- Here there are effectively two points to determine 
the topGOA fit. 

To test this approximation further, we have compared 
the topGOA quadrupole matrix elements with the matrix 
elements obtained by a full integration for a variety of 
nuclei and deformations. The result is shown in Fig. [31 
One can conclude that Eq. ^ is of sufficient accuracy 
for our purpose. 



D. Configuration mixing 



^{{q\RpQ2M)] 
= {q\Rp\q')t [ao + a2dlM + a, {d%{(i)f] , (9) 

where the coefficients of the polynomial depend on 
q, q' . As with the other matrix elements, it is important 
to include the point at (3 = 7r/2 when q and q' have 
opposite signs. 



As mentioned above, we typically compute about Nq — 
20 SCMF configurations to construct the energy land- 
scape. For many nuclei, however, only about half that 
number can be kept in the configuration mixing calcula- 
tion due to ill-conditioned norm matrices when the space 
is overcomplete. Nevertheless, the full configuration mix- 
ing calculation requires to compute about 50 projected 
matrix elements, which is beyond our computational re- 
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FIG. 2: The matrix element {q\Ri3Q2ii\<l) for the nucleus ^^Cr 
at an intrinsic deformation of q — 150 fm^. The open circles 
show the points used for evaluating the integral by a 12-point 
Gaussian quadrature. The solid circles show the points used 
for the topGOA, and the resulting fit. The three panels show 
the results for = 0, 1 and 2 going up from the bottom panel. 



20 
15 
10 
5 

-5 
-10 
-15 
-20 



0) 



' 1 

" □ 


1 1 1 1 1 

Q20 


' 1 ' 1 ' 1 ' y 


. o 


Q22 




- A 


Q22 




A 1 , 1 , 1 , 


,1.1,1. 



-20 -15 -10 -5 5 10 15 20 
GOA (c b) 

FIG. 3: Comparison of the topGOA and full calculation for 
Q21J. for various nuclei and deformations. 



culating quadrupole matrix elements between states of 
different angular momentum. The matrix element can 
change sign, depending on the deformations. For matrix 
elements connecting different manifolds of states (0+ and 
2+), there is no diagonal clement to anchor the GOA. 

Lacking a reliable GOA to determine the off-diagonal 
quadrupole matrix elements, we took another approach 
to simplify the configuration mixing calculation. The 
number of configurations has been reduced for each nu- 
cleus to a number small enough to make a global cal- 
culation feasible but large enough to have a sufficient 
accuracy on the energy of the lowest 0+ and 2+ states. 
Since we have to deal with energy curves of very different 
topologies, some care must be taken into the selection of 
points. The procedure that we have followed is explained 
in appendix [B] The number of selected configurations 
varies from 3 to 10, but is most often equal to 6. We 
have therefore labeled this approximation HW-6. 

The points selected for "^^Ar are shown in Fig. [T] as 
black circles. In this case, the single-configuration MAP 
energy of the ground state is E^olMAF) = -332.25 MeV, 
as quoted in the last section. The gain in energy from 
configuration mixing with the large set of configurations 
(11 in this case) is Eo{HW) - £;o(MAP) = -332.69 + 
332.25 = -0.44 MeV. This is to be compared with 
-332.61 -I- 332.25 = -0.36 MeV for the HW6 space. The 
error, 0.08 MeV, is within our targeted limit of accuracy. 
For the J — 2 projected wave functions, the energy gain 
by the HW treatment is —0.68 MeV and the error of the 
6-configuration truncation is 0.07 MeV with the same 
sign as in the ground state energy. With our present 
computer resources, we were able to test the HW6 trun- 
cation for about 100 nuclei. The 2+ excitation energies 
computed both ways are compared in Fig.|4l The approx- 
imation reproduces the energies to an r.m.s. accuracy of 
better than 0.1 MeV. The worst cases are two nuclei with 
coexisting minima at low excitation energy that are sep- 
arated by very low barriers, ^*^Pb and ^^''Pb, visible as 
points off the fine at the bottom left-hand corner of the 
graph. 

The same comparison for the reduced matrix element 
of the 2+ 0+ transition is shown in Fig. \5\ The agree- 
ment is very good except for the light Pb isotopes, ^^^Pb, 
^^^Pb, and ^^"Pb. Among the nuclides, the light Pb iso- 
topes are rather singular and we shall examine ^*^Pb in 
more detail in the next section. Overall, the accuracy of 
the HW-6 approximation is more than adequate for the 
present global survey. 



sources for a study of several hundred nuclei. In Ref. '23|, 
a GOA was developed for a coordinate corresponding to 
the deformation permitting calculations to be made to 
the needed accuracy only using the diagonal and subdi- 
agonal elements of the configuration mixing matrix, i.e. 
about 2Nq projected matrix elements. Unfortunately, the 
demands on the approximation are more severe when cal- 



III. SOME EXAMPLES 

In this section we shall examine the results for a sam- 
ple of nuclei with energy maps of different topologies: a 
light doubly-magic system, *°Ca; a heavy doubly-magic 
system, ^°^Pb; a transitional nucleus near magicity, '^^Ar; 
a soft nucleus exhibiting triple shape coexistence, ^^^Pb; 
and a well-deformed heavy nucleus, ^"'"Pu. We first ex- 
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FIG. 4: Comparison of 2+ excitation energies for the full and 
the HW-6 bases. 
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FIG. 5: Comparison of the reduced matrix element (0||Q2||2) 
for the 2"*" 0^ transition obtained within the full and the 
HW-6 bases. 



amine the results of the MAP approximation. 



A. MAP 

The MAP approximation is a variation after projection 
method within the very limited subspace defined by the 
axial quadrupole operator. For each J-value, one finds 
the configuration leading to the lowest energy, which thus 
could be different for J = and J — 2. The results for 
the observables of interest are presented on the first line 
of Table [ijtogether with the experimental data on the 
last line [2y|. In all cases, except ^*^Pb, the calculated 
2+ excitation energy is too high. For ^^Ar, this over- 
estimation has been attributed to the structure of the 
state @ . It is indeed very likely predominantly a broken- 
pair two-quasiparticle configuration, where the two occu- 
pied magnetic substates in the proton ^3/2 shell below the 
Z = 20 gap are coupled to J = 2. A rough estimate for 
its excitation energy is provided by two times the proton 



pairing gap, which leads to an excitation energy close to 
the experimental one. The description of this state re- 
quires the breaking of time-reversal reversal invariance 
and of axial symmetry in the SCMF, which is outside of 
what we can currently handle within our beyond-mean- 
field approaches. The excited state in the ''"Ca is of a 
different nature. Since, within the shell model, one does 
not obtain low- lying even-parity excitations, this state 
has been famous in the literature as an early example of 
shape coexistence. The first excited state of ^°Ca is a 0+ 
which is the head of an intruder deformed band. A de- 
tailed study of this nucleus with the full projected GCM 
has indeed obtained such a band [1] . Although the SLy6 
interaction was used in that case, the results that we find 
here are very similar. 

The next nucleus ^^^Pb shows still another kind of be- 
havior. It is a nucleus with several coexisting minima, 
which are separated by tiny barriers only. This isotope 
has been studied in detail in Ref. [7,] with several other 
neutron-deficient Pb isotopes. Under the MAP approxi- 
mation, the ground state is obtained from a configuration 
close to sphericity, while the minimum for the first 2+ is 
strongly oblate. The next example, ^°*Pb, is a heavy 
doubly magic nucleus. Simple shell model considerations 
give already an idea of what should be the dominant 
component of the first 2"*" excitation. It can be obtained 
by promoting a neutron from the occupied 113/2 shell to 
the unoccupied (79/2 shell, or a proton from the ft.11/2 
shell to the /ig/2 shell. The single-particle energy dif- 
ferences in the spherical mean-field configuration are 6.4 
and 5.9 MeV, respectively. The MAP energy is 2 MeV 
higher than the particle-hole energy. Again, like in the 
case of "^^Ar, the relevant configurations are broken-pair 
two-quasiparticle states outside our configuration space. 
The last nucleus in the table, ^^°Pu, is highly deformed. 
Its character is already seen in the SCMF wave func- 
tion, which has an intrinsic mass quadrupole moment of 
q = 3000 fm-^, of which 1145 fm^ is taken up by the 
electric quadrupole moment. Assuming that the wave 
function corresponds to a rigid rotor, Eq. (jAlip . one ob- 
tains a transition quadrupole moment (0||Q2||2) = 361 
e fm^ in agreement with experiment. The MAP approx- 
imation does not change matters; the minimizing q of 
the J — and J — 2 projected states are very close 
to that of the SCMF ground state. However, one sees 
from the table that the excitation energy of the 2+ state 
is too high by nearly a factor of two. This is another 
well-known problem, which has been seen in virtually all 
calculations using methods similar to ours: much bet- 
ter agreement would be obtained using the cranked HFB 
method to generate a wave function for the 2+ state. 



B. HW and HW-6 

We now examine the effects of configuration mixing on 
the properties of the 2+ state, which are also given in 
Table [D Mixing in the large ("full") configuration space 
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TABLE I: Results for some selected nuclei. The reduced tran- 
sition matrix element and the spectroscopic quadrupole mo- 
ment are defined in Eqs. IA2I and (|A8|) . respectively; the re- 
lation to the B{E2) is given in Eq. HA6|) . The theoretical 
values for ^***Pb from Ref. assume that the ground state is 
spherical and the excited state is oblate. 
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0.083 


377 


-340 
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significantly reduces the energy in two cases, "Ca and 
^°®Pb, raises the energy in one case, ^*^Pb, and has httle 
or no effect in two cases, ■^^Ar and ^^°Pu. The insensitiv- 
ity for ^^°Pu is to be expected since strongly deformed 
rotors do not have large shape fluctuations, see the de- 
tailed discussion for this example in Ref. Q. Including 
shape fluctuations improves the description of a soft nu- 
cleus as i®*Pb. We see that the full calculation (HW) 
produces an excitation energy that approaches the ex- 
perimental value. In all the cases where the fluctuations 
change the energy, the change goes in the right direction 
and decreases the theoretical error. 

On the third line of Table U we show the effect of the 
truncation of the configuration space in the IIW-6 ap- 
proximation. In all but the case of ^^*Pb the energies 
are close to the full HW results. The light Pb isotopes 
are quite exceptional, but we saw in the previous sec- 
tion that HW-6 is reliable enough for a global survey. 
The next line in Table U shows results from other calcu- 
lations. The reported calculations of ^^^Ar, *°Ca,^^^Pb 



and ^ °Pu were done with the full projected GCM with- 
out approximations, using the same computer codes as 
here, but with a slightly different energy functional. We 
see that the results are qualitatively similar to what we 
found, indicating a mild sensitivity to the specific energy 
functional. 

We now discuss the quadrupole matrix elements in 
more detail. The simplest case is ^^°Pu, which, as dis- 
cussed above, behaves very much like a rigid rotor. In 
the rotor limit, the transition quadrupole matrix element 
is proportional to the spectroscopic quadrupole moment 
of the 2+ state. The relation is given in appendix A, 
(3c/(0||(32||2) « 0.9. We see from the table that this is 
well satisfied for all calculations of ^'^''Pu. For the non- 
deformed nuclei, the spectroscopic quadrupole moment 
is small, as would be expected for a spherical vibrator. 
In the four cases given in Table [H the HW and HW-6 
transition matrix elements, although overestimating the 
experimental data, are better than a factor two, even in 
cases where the dominant component of the 2^ appears 
to be incorrect. This is probably related to the fact that 
the quadrupole moment is a bulk property that is en- 
tirely determined by the overall distribution of the local 
density, while the energy is sensitive to the detailed struc- 
ture and occupation of each single-particle wave function. 
Note also that allowing spreading of the wave functions 
over several configurations improves the MAP result. 



IV. GLOBAL PERFORMANCE 

We carried out the MAP and HW-6 calculations on 
even-even nuclei with known binding energies, excluding 
light nuclei with iV or Z < 8. This is the set studied in 
Ref. Of these, 522 have known 2+ excitation ener- 
gies. These energies range from 39 keV to 6.9 MeV, thus 
spanning more than 2 orders of magnitude. The theo- 
retical numbers span the same range, but as we saw in 
the last section there can easily be a factor two error in 
specific cases. 

In view of the results of the previous section, we have 
excluded from the full set of nuclei the ones for which 
one can have suspicion about our approximation scheme. 
To identify these nuclei, we have compared our present 
HW-6 results with the global calculation performed ear- 
lier where the number of configurations included in the 
calculation of the ground state was not limited. We elim- 
inate all the nuclei for which the difference between both 
calculations for the energy of the 0+ ground state was 
larger than 250 keV. The selected set of nuclei does not 
include ^^^Pb, and similar nuclei which are too soft to be 
represented by either a MAP calculation or a small num- 
ber of quadrupole configurations. Out of the 522 nuclei 
calculated, there remain after selection 359 . 
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TABLE II: Statistics for the performance of the MAP and 
HW-6 calculations. 
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FIG. 6: Histogram of the logarithmic errors Re of the energies 
of the first excited 2+ states (solid) and the logarithmic errors 
Rq of the reduced matrix elements (0||(32||2) (dashed) for the 
359 even-even nuclei included in our survey (solid). The two 
panels show the results of the MAP theory (left) and the 
projected GCM in the HW-6 approximation (right). 



A. Global results 

Because the energies span a broad range and the error 
can be large, we will quote the aggregated results for the 
logarithm of the ratio of the theoretical to experimental 
energies, 

Re = \og{Eth/E,^p). (10) 

A histogram of this quantity for the entire set of nuclei 
is shown in Fig. [51 displaying the MAP results in the 
lefthand panel and the HW-6 results in the righthand 
panel. We see that the results of both methods tend to 
be too high, with a fairly broad distribution containing 
both negative and positive errors. Quantitative statisti- 
cal measures of the distribution are given in Table HIl 

The average MAP error is found to be (Re) ~ 0.28 
but the average of the absolute value of the error is much 
larger, {\Re\) = 0.48, corresponding to an error of the 
order of 66 %. The dispersion around the average is 
also quite large: {{Re - {R))Y^^ = 0.49. With such 
a dispersion, an error larger than a factor of two is not 
unusual. Specifically, of the 359 nuclei in the data set, 
19 % have a calculated energy too large by a factor of 
two and 4 % are too low by the same factor. Fig. [7] 
shows a scatter plot of the MAP and the experimental 
energies. One sees that the energies are overestimated for 
most nuclei, and in particular for nuclei with either a low 
or a high excitation energy of the 2"*", where the nuclei 
are predominantly deformed or magic, respectively. This 
is consistent with what we saw in the examples of the 
previous section. For excitation energies in the range 
200 keV to 1 MeV, there is no obvious trend in error of 
the MAP calculation. 

As can be seen in Table |TT1 the mean error of the HW- 
6 calculation is significantly larger than the MAP error, 
with an average around {Re) ~ 0.51. The dispersion 
around the average is however lower and the average of 



the absolute value of the error is only slightly larger than 
for the MAP results (0.54 compared to 0.48). In Fig.[7]are 
plotted the HW-6 excitation energies as a function of the 
experimental data. One can see that in most cases, the 
2+ excitation is overestimated and this tendency is much 
more pronounced than for the MAP results. In fact, in 
many cases, the HW-6 2+ energy is larger than the MAP 
one. This increase when the configuration mixing corre- 
lations are included means that the correlation energies 
predicted by our method are larger in the ground state 
than in the 2+ state. There can be many origins to this 
difference in correlation energies. The lack of triaxial con- 
figurations certainly affects more deeply the states with 
J ^ since for these states the spherical point does not 
contribute and the coupling between prolate and oblate 
configurations is disfavored. It is also clear that the MAP 
procedure is better defined numerically than the HW-6 
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FIG. 7: Scatter plot comparing the theoretical and experi- 
mental 2^ excitation energies of the 359 nuclei included in the 
survey. The two panels show the results for the MAP method 
on the left and the HW-6 approximation on the right. 
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one. In each case, we are sure to have determined for 
both J = and J — 2 the quadrupole moment giving 
the minimal energy after projection. For the configura- 
tion mixing, the fact that we have excluded the nuclei for 
which the J = energy is too different from our previous 
global calculation makes the determination of the 0+ en- 
ergy reliable. We do not have a similar check for J — 2 
and there are cases where the number and the spacing 
of points taken for J = 2 is not fully adequate and the 
energy of this state less accurate. 

While the energies are not accurately predicted, the 
quadrupole properties come out much better and with 
rather similar errors for both the MAP and HW-6 results. 
It is well known that the intrinsic quadrupole moments 
of deformed nuclei are rather insensitive to the details 
of the energy functional, and indeed we found that the 
quadrupole transition matrix element is much better de- 
termined overall than the energy. The dashed histograms 
in Fig. [6] show the logarithmic ratio Rq of the reduced 
quadrupole transition matrix elements (0||(52||2) (see ap- 
pendix The average error is only 0.12 for the MAP 
calculation, corresponding to matrix elements that are 
15% too large and 0.09 for the HW-6 results (error of 
9%). The r.m.s. spread is also reduced. For example, in 
the MAP case, it takes the value 0.22 corresponding to 
transition matrix elements that are -10% to -1-46% of the 
data. We will now analyze separately the different kinds 
of nuclei. To that aim, we divide the nuclei by type and 
examine in more detail the performance for subgroups 
that are deformed, doubly-magic, and singly-magic. 
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FIG. 8: 2"*" excitation energies in deformed nuclei as a function 
of neutron number: MAP in left panel; HW-6 approximation 
in right panel. Experimental data are shown as diamonds. 
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1. Deformed nuclei 

We first have to define a criterion to select which nu- 
clei should be considered deformed. Obviously, there is 
no rigorous division of nuclear types, and any division is 
somewhat arbitrary. One possibility is to make a selec- 
tion on the basis of the intrinsic quadrupole moment of 
the MAP ground state, taking into account overall size 
effects by using the geometric shape parameter P2 [de- 
fined in Eq. (|A13p ] to make the selection. This criterion 
will catch many light nuclei along with the usual nuclei 
in the lanthanide and actinide regions. One should add 
the criterion of rigidity to the selection as well to elimi- 
nate the nuclei that have large fluctuations in shape. In 
this sense, what we are seeking to categorize are nuclei 
that behave like rigid rotors. A criterion that makes a 
nice selection is to demand that the average deformation 
(32 is large than the r.m.s. fluctuation about the average, 
P2 > ((/32 — (32)^^^- These quantities are computed using 
the full HW wave functions of Ref. [3| , and the criterion 
selects 134 deformed nuclei from our set of 359. Their 
energies are plotted as a function of neutron number in 
Figs. [S] with the MAP results in the lefthand panel and 
the HW-6 results in the righthand panel. The two plots 
are rather similar. 

We see that the predictions are too high for the ac- 



FIG. 9: Ratio of theoretical and experimental transition 
quadrupole matrix elements (0||(52||2) in deformed nuclei, as 
a function of neutron number. 



tinides, while on the average they are quite reasonable 
for rare earth nuclei. The statistic on the errors for de- 
formed nuclei is summarized in Table [ill One can see 
that the average error is smaller than for the full set. 
The dispersion in the error is the same for both the MAP 
and the HW6 approximations, so, the axial quadrupole 
correlations do not seem to be the source of the nucleus- 
to-nucleus fluctuations of error. 

Fig. [9] shows the ratio of theoretical to experimental 
quadrupole transition matrix elements for the deformed 
nuclei. Here the actinide nuclei come out very well. 
There is more fluctuation in the rare earth and the light 
nuclei that qualify as strongly deformed but the overall 
results are quite satisfactory. 



2. Magic and semi-magic nuclei 

We now turn to doubly- and singly-magic nuclei, which 
present quite different problems for the theory. The com- 
parison between theoretical and experimental 2+ excita- 
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FIG. 10: 2^ excitation energies near magic nuclei, as a function of neutron and proton numbers. MAP and HW6 results are 
shown by points connected with lines, while experimental data are shown by filled diamonds. 
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FIG. 11: Ratio of theoretical and experimental quadrupole 
transition matrix elements near magic nuclei as a function 
of neutron and proton numbers. MAP and HW6 results are 
shown as points connected by dotted and solid lines, respec- 
tively. 



TABLE III: Excitation energy of the first 2+ state in MeV for 
doubly magic nuclei. 



N Z exp. MAP GCM 



20 20 3.9 

28 20 3.8 

28 28 2.7 

82 50 4.0 

126 82 4.1 



5.4 
2.7 
2.0 
6.3 
7.9 



5.4 
2.7 
2.2 
5.8 
6.7 



tion energies of six doubly magic nuclei is shown in Ta- 
ble mil The MAP and HW-6 energies are too high in 
three cases and too low for the other three, preventing 
us from drawing any general conclusions. 

There are 71 semi-magic nuclei in our compilation, re- 
quiring that either the neutron or the proton number 
equal to 28, 50, 82 or 126. Graphs of 2+ excitation en- 
ergies are shown in Figs. [TO] as a function of neutron 
number and proton number. MAP results are shown 



by the points connected with solid lines, HW-6 results 
with long dashed lines, and experiment by shorter dashed 
lines. For both MAP and HW-6, the excitation energy 
has a peak at the doubly-magic nuclides which decreases 
gradually going away from that nucleus. In contrast, the 
experimental peak is a sharp spike at the doubly-magic 
nuclides. This deficiency of the theory is related to the 
absence of broken-pair two-quasiparticle excitations that 
can give a lower non-collective 2+ state as we saw already 
in the example of '^'^Ar. As to be expected from this dis- 
cussion, the statistical measures are much poorer for this 
class of nuclei. The average calculated energy is 50% 
higher than experiment and the average calculated tran- 
sition quadrupole moment is 2.2 times the experimental 
value. 

Fig. [TT] compares the theoretical and experimental 
transition quadrupole moments in scmimagic and magic 
nuclei. The data is much more meager than for the ener- 
gies, but one can see that the theory nearly always is too 
high. As discussed earlier, this is to be expected when 
the lowest 2+ is not collective. 



B. Discussion 

In view of their restricted form and of the way they 
have been fitted to selected experimental data, the cur- 
rent energy functionals are certainly too limited and, as 
discussed in [1, |2^ , present deficiencies which are at the 
origin of some of the discrepancies between our calcu- 
lations and experiment. However, the present analysis 
clearly points also to deficiencies of the variational space 
that is used, that affect more excitation energies than 
quadrupole moments. 

One can expect that our configuration space spanned 
by axial quadrupole SCMF wave functions covers the cor- 
relations that dominate the description of the 0"*" ground 
states of even-even nuclei. On the contrary, there are sev- 
eral competing possibilities to construct a low-lying 2^^ 
state, some of them being completely absent from our 
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description: 

1. a broken-pair two-quasiparticle excitation within a 
partly-filled j shell with an excitation energy of 
about two times the pairing gap (for example near- 
magic nuclei like ^^Ar) 

2. a broken pair two-quasiparticle excitation involving 
two different j shells, one occupied, the other unoc- 
cupied, with an excitation energy of about the gap 
between the j shells involved (for example doubly- 
magic nuclei as ^°*Pb) 

3. a collective vibrational state 

4. a collective rotational state (for example for well- 
deformed rare-earth and actinide nuclei) 

The states corresponding to such pure configurations 
should of course be mixed in actual nuclei. The projected 
GCM of axially deformed SCMF states that correspond 
to HFB vacuua, as it is used here, cover only the latter 
two of these configurations. In nuclei where the lowest 
2"*" state is dominated by broken-pair two-quasiparticle 
states, the 2+ state that our method enables us to de- 
scribe corresponds to a higher-lying collective 2"'' state. 
On the other hand, the number of nuclei for which the 
lowest 2+ state is indeed dominated by a broken-pair two- 
quasiparticle state can be expected to be small, and re- 
stricted to the immediate vicinity of doubly-magic nuclei. 
In nuclei where the first 2+ is expected to be collective, 
either vibrational, or rotational, the excitation energies 
are also on the average too high. This result confirms 
on a large scale previous studies performed with similar 
methods for smaller sets of nuclei [1, 0, B i, E lU [11 

An obviously missing degree of freedom is triaxiality, 
as our configuration mixing contains only two out of the 
five degrees of freedom of the quadrupole tensor. On 
the one hand, it is well-known that many transitional 
nuclei are 7 soft. A recent global study based on a semi- 
microscopic method has even indicated that the potential 
energy surface of many transitional nuclei might have a 
triaxial minimum f27l| , although with an energy gain that 
remains very small. It has also been shown [28| that 
in some cases triaxial quadrupole configurations can be 
more favorable after angular-momentum projection than 
axial configurations. The effect of triaxiality on excita- 
tion energies has also been studied with the help of an 
effective five-dimensional Bohr Hamiltonian derived from 
mean- field calculations using the Gogny force [l^,!!^ or a 
Skyrme interaction [3l| . A similar overestimation of the 
lowest 2"*" energy as in the present study has been found. 
Therefore, if the effects of triaxiality are certainly non- 
negligible on total energies, it is unclear whether they 
will improve excitation energies. 

The effect of triaxial quadrupole deformations on exci- 
tation energies is not obvious. The method that we use 
has however clearly an artifact which favors the 0"*" en- 
ergy with respect to the 2"*": the mean-field is optimized 
to describe the ground state and not excited states. An 



TABLE IV: Excitation energy of the first 2+ state in MeV for 
Zn isotopes. 



N Expt. MAP HW-6 cranked SCMF 



30 


0.89 


1.85 


1.89 


1.35 


32 
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34 
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36 
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0.37 


38 
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1.88 


1.95 


0.41 


40 


1.50 


2.64 


2.33 


0.80 



obvious improvement would be to perform an exact vari- 
ation after projection (VAP), performed separately for 
J — and J = 2. Variation after projection on angular 
momentum starting from effective interactions and with 
a full model space does not seem however within com- 
putational possibilities in a near future. A more mod- 
est approach would be to use the self-consistent cranking 
method to optimize separately the intrinsic configura- 
tions describing different spin states. This can be done 
by introducing in the mean-field equations a constraint 
on the projection of the angular momentum. It has been 
shown that the self-consistent cranking method is an ap- 
proximation of a variation after projection on angular 
momentum for well deformed nuclei [s^ . Such an ap- 
proach will necessarily improve the description of the 2+ 
states and decrease their excitation energy. 

Even for nuclei which are not deformed, the effect of 
a cranking constraint, which breaks time-reversal invari- 
ance can only go in the right direction, although it is not 
clear a priori that the introduction of two-quasiparticle 
excitations would not be even more important. To give 
some insights into the effect of a cranking constraint for 
isotopes which are not deformed, we show in Table HVl the 
results of MAP, HW-6 and cranked SCMF calculations 
together with the experimental data for some neutron- 
rich Zn isotopes. The cranked SCMF calculations were 
done with the method described in Ref. [S^ and with the 
same effective interaction as for the other calculations of 
this study. Both MAP and HW-6 results strongly over- 
estimate the experimental energies, the MAP being even 
rather irregular in its predictions. The cranking results 
are given in the last column. The 2^^ excitation energy in 
this case is the difference between the mean-field ground 
state energy and that of the state obtained with a crank- 
ing constraint = 2. One sees that the energies are 
significantly lower compared to the calculations where 
time reversal invariance is imposed. The numbers that 
are obtained can even be lower than the experimental 
data. One cannot go too far in the interpretation of these 
results which do not include any projections. However, 
they show that an optimization of the 2"*" wave function 
with a cranking constraint might have a significant effect 
going in the right direction for all nuclei. 

It remains to verify what will happen when exact pro- 
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FIG. 12: Chart of nuclides showing the even-even nuclei for 
which the excitation energy of the first 2^ state is known. 
Nuclei for which the MAP calculation is in error by more 
than a factor of two are shown by solid triangles. Upper 
panel; theory too low; lower panel: theory too high. 



FIG. 13: Chart of nuclides showing the even-even nuclei for 
which the transition quadrupole moment to the first 2"*" state 
is known. Nuclei for which the MAP calculation is in error by 
more than a factor of two are shown by solid triangles. Upper 
panel: theory too low; lower panel: theory too high. 



jection and configuration mixing will be performed, but 
the perspective of a global qualitative improvement of 
the present results seems reasonable. The generaliza- 
tion of the method used here for use with cranked SCMF 
states requires a similar, even greater, effort as the gen- 
eralization to triaxial shapes: the cranking constraint in- 
duces non- axial intrinsic currents, even when the over- 
all density distribution remains axial. The broken-pair 
two-quasiparticle states discussed above require a similar 
generalization of the spatial symmetries in the projected 
GCM as the cranked HFB vacuum. 



V. SUMMARY AND OUTLOOK 

This is the first systematic survey of the 2+ excita- 
tions of even-even nuclei in the framework of a mean- 
field based method including symmetry restoration and 
starting from an effective energy functional. The effort 
necessary for this task is significantly larger than the one 
required for our earlier study on the ground states of 
these nuclei ^ , both for the representation of the collec- 
tive 2"*" states, and to obtain a sufRcient precision for the 
matrix elements of the quadrupole operator. For several 
nuclei we could not reach an acceptable precision, so that 
the subsequent analysis was performed on a reduced set 
of data. 

Qualitatively the excitation energies and B{E2) val- 



ues track the data for the great majority of the 359 nu- 
clei studied. However, predicted energies and B{E2) val- 
ues are systematically too high, and there are a number 
of cases where key ingredients are clearly missing. The 
worst cases, where the observable is more than a factor 
of two in error, are marked on the charts of nuclides in 
Figs. [13 and [T31 One can discern some patterns that 
point to deficiencies in the energy functional and in the 
GCM methodology that both may be correctable. 

Many of the outlying points in Fig. [12] are cases where 
the theory predicts a nearly spherical nucleus while the 
data show it to be deformed, or vice versa. An example is 
*'^Zr, predicted spherical but obviously deformed in view 
of the very low excitation energy of its first 2+. The 
shell effect predicted by our effective interaction for ^°Zr 
is clearly too large, as already analyzed in our study of 
ground-state correlations Q . All conclusions of Ref . [3| 
about necessary future work on the effective interactions 
and the model space also apply here, see also [25j . 

Our calculation reproduces rather nicely the 
quadrupole transition matrix elements between the 
first 2+ and the ground state: the average error that we 
obtain is around 25%. The situation is less satisfactory 
for the excitation energies of the first 2+ states, which 
is nearly always overestimated. This seems to be a 
general problem that has been noticed before in many 
calculations using Skyrme and Gogny interactions. As 
argued above, we relate this deficiency mainly to the 
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current restrictions of the variational space that we 
use. To overcome this hmitation, the extension of the 
variational space to include triaxial states, and cranked 
SCMF states is highly desirable. Work in that direction 
is underway. The enormous increase in computational 
time, however, will not permit its large-scale application 
right away. 
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APPENDIX A: MATRIX ELEMENTS OF THE 
QUADRUPOLE OPERATOR 

For reference, we quote the definitions of the 
quadrupole matrix elements and simplified versions of 
the formulas from Ref. 6] for calculating them. For 
the sake of simple notation, we will give all expressions, 
where applicable, for matrix elements between two dif- 



ferent SCMF configurations, \q) and \q'). The general- 
ization to GCM states with their weighted summation is 
straightforward as it does not affect the angular momen- 
tum algebra. 

The electric quadrupole operator is defined as 

Q2^^eJ2r^Y,^ih) (Al) 
p 

We start with Eq. (A7) of Ref. for the reduced ma- 
trix element of the quadrupole operator between two pro- 
jected axial states 

{MQ2\\JW) = ^ E {J'02^^\J^,) 

X / rfcos(/3X(/3) (<i\RfiQ2M) 
Jo 

(A2) 

where J and J' are assumed to be integer and even. Rp 
is the rotation operator, d'l^^{(3) is the Wigner c?-function, 
and Afjq is the normalization of the J-projected SCMF 
state, Eq. The reduced matrix element on the left- 
hand side is defined as [s^l 

{JMq\Q,,\J'M'q') ^ il^I^^{Jq\\Q,\\j'q'). 

\' ZJ + i 

(A3) 

Equation (jA2| can be even simplified further using the 
symmetries of the Wigner functions and the quadrupole 
operators 



{Jq\\Q2\\J'q') = 



V2J+l{2J' + 1) 



J'q' 



J'020|J0) / dcos(/3) d^oW {q\Ri3Q2o\q') 



+2j2{J'02fi\Jfi) / dcos(/3) d'o^iP) 5R{(g|i?/3Q2Ml'z')} , 



r 



(A4) 



which serves as the starting point for the GOA set-up in 
section III CI 

To compute the matrix element for the 2+ — > 0+ tran- 
sition, one evaluates the above formula with J = and 
J' = 2. Only fi = contributes in this case and the result 



matrix element by [34 



B{E2,2^ 



(0g||Q2||2g')^ (A6) 



The other matrix element of interest is the spectroscopic 
quadrupole moment of the J = 2 excited state, defined 
as 



(0g||Q2||2g') 



V5 



A/'2gA/'og' Jo 



dcos(/3) {q\RpQ2o\q'). 



(A5) 



Q, = ^l^{22q\Q2o\22q') 



(AT) 



The B{E2) for the transition is related to the reduced In this case the sum over cannot be avoided. The final 
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result is 



Qc 



Mi 



-iMo + ^Ah + ^Ah] (A8) 



where 



16n 



dcos{p)dl^{l3){q\RpQ2M')- (A9) 



The rotor model provides a convenient reference for 
estimating quadrupole matrix elements. In terms of 
the intrinsic quadrupole moment of the configuration, 
(g I (320 1?), the relations are 



Qc 



167r 



q\Q2o\q) 



and 



(0g||Q2||2<z)rotor = {q\Q2a\q) 



(AlO) 



(All) 



Finally, we specify the deformation of a configuration by 
the mass quadrupole moment with the spectroscopic nor- 
malization. The relation is 

q = \{q\Y.rnico^He)-l]\q). (A12) 



We also use the dimensionless deformation parameter /32 
defined by the equation 



/32 = 



(A13) 



3 ARl' 

using the liquid drop radius constant Rq = 1.2 A^^^ fm. 



APPENDIX B: SELECTION OF 
CONFIGURATIONS IN HW-6 

In this appendix we describe in more detail how the 
configuration set was chosen for the configuration mixing. 



The following rules were applied to select configurations 
for each nucleus and for J = and J = 2. The rules are: 



• Start with the set of 15-20 constrained configura- 
tions that were used in our previous study [J]. 



• Divide the set into prolate and oblate configura- 
tions. For both sets and each angular momen- 
tum value, find the configurations |(?rnin) that have 
the minimum energy after particle-number and 
angular-momentum projection. 



• In each set, select the projected configurations on 
each side of the minima that have an overlap close 
to, but larger than 0.5 with |gmin). This leaves both 
sets with up to three configurations. 

• Join the prolate and oblate sets, taking out the 
oblate configuration with the lowest deformation 
if its overlap with the least deformed prolate con- 
figuration is greater than 0.9. 



• Add to the set of J = configurations all the J — 2 
configurations that do not overlap a J = config- 
uration by more than 0.9. Likewise add J ~ 
configurations to the J = 2 set. 

Most resulting sets include 5 to 6 configurations, al- 
though some could be larger or smaller. For example, for 
nuclei in the rare-earth and actinide regions that present 
a deep and narrow prolate minimum in the total energy 
surface, the oblate configurations are too high in energy 
to play a role and 3 configurations are sufficient. For 
some other nuclei, several points are needed to connect 
the prolate and oblate sets, making the configuration set 
larger than 6. In a few cases, the selected sets lead to 
instabilities in the solution of the HW equations, related 
to too small eigenvalues of the norm kernel. These cases 
had to be treated by hand to select the configurations. 
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